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Abstract 

A new method ( PI-DFT ) which combines path integrals and density func- 
tional theory is proposed as a pathway to many fields of physics. Within path 
integral theory it is possible to construct particle densities without explicitly 
calculating individual wave functions. These densities can directly be used as 
an input to energy density functionals. Thus our method makes full use of 
the theorem of Hohenberg, Kohn and Sham which shows, that the energy of 
a many electron system only depends on the particle density. 
At glance we present a recursion formula for the calculation of many fermion 
and boson particle densities from one-particle densities at a set of different 
temperatures. For both statistics the numerical effort of our method increases 
only with the square of the particle number. 
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Path integrals and Density Functional Theory have both been investigated and developed 
extensively in the past decades [|. Here we present a new method (PI-DFT), which formally 
joins both fields and for the first time uses the theorem of Hohenberg, Kohn and Sham |^,^ 
to its full extent. While the method is applicable to bosons as well, the most interesting case 
is that of fermions because of its great relevance in quantum chemistry. The pride of our 
method is that we never deal with individual wave functions to compute the particle densi- 
ties but still solve the full problem. The basic underlying idea is a self-consistent iteration 
scheme, in which the particle density is computed by means of path integrals, while effective 
potentials and observables are computed with methods lent from density functional theory. 
For the problem of constructing completely antisymmetric or symmetric particle densities 
we propose a recursion formula which has originally been invented by us for the partition 
function [^]. With this formula there is no need to know the individual wave functions 
to construct the appropiate particle densities for Fermi-Dirac (FD) or Bose-Einstein (BE) 
statistics. Instead, these are calculated from single-particle densities but at different tem- 
peratures. 

Using this formula the numerical effort for both statistics grows only with the square of the 
particle number. 

In the following we will consider N-particle hamiltonians of the from 

i.2 N j2 N N 

2m ^ df2 ^ 

That is, we assume the potential to contain one and two particle interaction terms as en- 
countered typically in quantum chemical problems. 

In conventional Path Integral Monte Carlo techniques for fermion systems an NM di- 
mensional integral with an integrand containing an ( x A^) determinant has to be calculated 
( see e.g. resulting in a numerical effort proportional to A^'^. Equation (|^) can be de- 

^See for example [Q,^ and references therein. 



2 



rived easily using the Trotter formula [§] and is equivalent to the thermal Path Integral as 
the number of time slices M goes to infinity. 
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Although the method is suitable for systems like ^He clusters and most physical relevant 
observables can be calculated, its application is numerically very extensive in general. 



According to the theorem of Hohenberg, Kohn and Sham and its extension to 
thermal systems given by Mermin the potential given by the term in square brackets in 
(0) can be replaced by an effective one-particle potential. 

(p{x) = v{x) + / dx'ri{x') u{\ x' — X \) . (3) 



Since the energy is only a functional of the one particle density r]{x), which can be easily 
calculated within the approximation (|^) by Monte Carlo methods and f{x'), the solution 
of the stated problem, can be found by self-consistent iteration of (H) and (||). A priori the 
use of this approximation has no advantages compared to the direct calculation of (H), but 
together with a much more effective way of calculating the partition function, which will be 
stated in the following, the above iteration scheme develops to a very powerful tool. 

Let us denote with r]N{x; (3) the probability of finding a particle at position x in a 
system of N noninteracting fermions or bosons moving in a common potential thermalized 
at temperature (3 and with Z^^jS) the corresponding canonical partition function. For such a 
system it is sufficient to compute the one-particle densities ?7i(x; ?7i(af; 2/3), . . . , 77i(x; 
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but at a set of different temperatures. Then Zn{(3) and r]]^{x,[3) can be calculated by the 
recursion formulas 

1 ^ 

Zn{(3) = ^ Y.i±)''''Zim Z^_,{(3) (4) 
fc=i 

N 

Vn{x; P) = Z],\P) J2i±)'+'Z^ik(3) 771 (f, kp) Z^^M (5) 

k=l 

where the plus and the minus signs stand for Fermi-Dirac and Bose-Einstein statistics, 
respectively, and Zo{P) = 1. 

The proof for @ has already been given explicitly in The validity of (^) can be proven 
in the same manner. As the whole procedure is somewhat lengthy, we omit it here and 
sketch the basic idea using the two-fermion system as an example. 

Let efc denote the energy eigenvalues of the one-particle system and ipk the corresponding 
eigenfunctions. Then the partition function is given by 

Z2{f3) = WY.^^Pi-f^i^>^ + ^i)) ■ (6) 

^ k lytk 

The probability pk of finding a particle in the one-particle state | k) is simply 

p, = Z^\(3) I ^exp(-/?(6, + ei)) (7) 

= Z2\P) I Eexp (-/5(efc + q)) - exp {-2/3ek) 
= Z^\(3) i [exp(-/5e,) Z^{(3) - exp(-2/?efc)] 
For the particle density this yields 

r]2{x](3) = J2Pk I i^k{x) P (8) 

k 

= Z^\f]) I Mx) r [exp(-/?efc) Zi(/9) -exp(-2/5Q)] 

k 

= Z^\P) [z^iP) ?7i(f;/3) - Z,{2P) vi{S;2P)] , 
which is nothing but equation for the special case N = 2. 

A simple example shows that the given formulas work very well and exactly. We calculated 
the exact solvable system of uncoupled harmonic oscillators. 
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for both statistics. The one particle densities kl3) shown in Fig. 1 are calculated with 
standard path integral techniques (see e.g. In Fig- 2 and Fig. 3 the exact particle 

densities and the densities calculated via are shown for both statistics. The outcome is 
an almost perfect agreement with the exact results. The computed energies agree equally 
well with the exact values. For our examples the errors were in all cases below 0.1 %. 
Obviously it is essential in the fermion case, that there is a nonvanishing occupation prob- 
ability of the excited states for the chosen temperature, i.e. the results of the one-particle 
Path Integral Monte Carlo simulations have to contain information about the excited states. 
Because of the fact that in the low temperature range the energy is in general a very slowly 
increasing function in the case of Fermi-Dirac statistics, this is not a too bad restriction, 
even if one is going to calculate ground state properties of the fermion system. 

Possibly the most interesting topic concerning PI-DFT is an analysis of its computational 
costs. The basic steps in PI-DFT are the following: 

1. Construction of an initial guess for the N particle density. 

2. Calculation of the effective potential. 

3. Calculation of single-particle densities and the corresponding energies at N different 
temperatures. 

4. Calculation of the partition function by spline interpolation of the caloric curve along 
with numerical integration using 
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Calculation of the N-particle partition function using (^. 
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Calculation of the N-particle density using (|]). 
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Steps 2-6 have to be repeated until the particle densities and the energy are converged. The 
main numerical effort lies in the steps 2 and 3. The calculation of the effective potential 
is usually proportional to iV^. Because of the fact that the numerical effort of a Path 
Integral Monte Carlo simulation grows normally proportionally to (5 the same applies to 
step 3. Steps 4 and 5 are negligible, while the effort for step 6 again grows proportionally 
to N"^ assuming that the number of used grid points is approximately proportional to the 
number of particles. Overall this indeed yields a numerical effort proportional to the square 
of the particle number for PI-DFT. Hence it follows that PI-DFT can be seen as a formal 
proof, that e.g. quantum chemical calculations can be done with such a slowly increasing 
numerical effort. We hope that this will encourage further improvements of other related 
methods, too. This paper presents only the onset of PI-DFT and a lot of improvements seem 
to be possible and required. For example, the recursion for the partition function might be 
numerically unstable under unfavorable circumstances in the case of FD statistics because 
of the alternating sums. 

In a forthcoming paper we will present a study on '^He clusters as a first physical application 
of PI-DFT. 
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FIGURES 

FIG. 1. Particle densities r]\{x\f3) calculated with Path Integral Monte Carlo for various tem- 
peratures (Ji = k = m = 1). 

FIG. 2. Exact harmonic oscillator particle densities and particle densities calculated with Path 
Integral Monte Carlo in connection with equation (3) for N=2 at T=0.5 and N=4 at T=2.0 for 
the case of FD statistics {Ti = k = m = l). 

FIG. 3. Exact harmonic oscillator particle densities and particle densities calculated with Path 
Integral Mont Carlo in connection with equation (3) for n=2 at T=0.5 and N=4 at T=1.0 for the 
case of BE statistics {li = k = m = l). 



8 



T=2.0 

T=1.0 




— N=4, T=2.0 ( exact) 
-- N=4, T=2.0 ( PI-DFT) 

— N=2, T=0.5 (exact) 




10.0 



N=4, T=1 .0 (exact) 
N=4, T=1 .0 (PI-DFT) 
N=2, T=0.5 (exact) 
N=2, T=0.5 (PI-DFT) 
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